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Abstract 

We have observed the rare decay B° — > p + p~ in a sample of 89 million BB pairs recorded with 
the BABAR detector. We measure the branching fraction B(B° —> p + p~) = (2TZQ-r) x ^ & 
and determine the longitudinal polarization fraction Ti/T = 0.99^QQy ± 0.03. Our results are 
preliminary. 
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1 INTRODUCTION 



Charmless B meson decays provide an opportunity to measure the weak-interaction phases arising 
from the elements of the Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing matrix 1 . There 
has been increasing interest in the study of B — > -kit and pir decays where the time-dependent 
CP-violating asymmetries are related to CKM angle a = arg [ — V^VZ / V ud y* b ] and interference 
between the tree and penguin amplitudes could give rise to direct-CP violation. The decay 7 B° — > 
p + p~ is another promising mode for CP violation studies and constraints on the angle a. The 
angular correlations in this decay to two vector particles introduce additional complications in the 
analysis, but the measurement of the magnitudes or phases of the helicity amplitudes will provide 
better understanding of the decay models [2|Ei|[3|- 

The decay B° — > p + p~ is expected to proceed through tree-level b —* u transitions and CKM- 
suppressed b — ► d penguins as illustrated in Fig. The presence of penguins and both CP-even 

W+ 




Figure 1: Two of the diagrams describing the decays B° — > p + p . 



(S- and D-wave) and CP-odd (P-wave) components in the decay amplitude complicate the mea- 
surement of a. Isospin relations among the three B — » pp modes may reduce the uncertainties in 
the measurement of a due to penguin contributions (penguin pollution), analogous to the methods 
proposed for B — > tttt [3]. The recent limit on the B° — ► p°p° decay rate jH] and the measurements 
of the B + — > p + p° branching fraction Ej result in experimental limits on the amount of penguin 
pollution, while B + — ► polarization measurements provide evidence that the CP-even longi- 
tudinal component dominates in the B — > pp decay amplitudes. In this paper we report on the 
observation of the B° — ► p + p~ decay mode and the measurement of its branching fraction and the 
longitudinal polarization fraction in the decay. 

2 THE BABAR DETECTOR AND DATASET 

In this analysis, we use the data collected with the BABAR detector [S] at the PEP-II asymmetric- 
energy e + e~ collider |Hj operated at the center-of-mass (CM) energy of the T(4S) resonance (y/s = 
10.58 GeV). These data represent an integrated luminosity of 81.9 fb _1 , corresponding to 88.9 
million BB pairs, at the T(4S) energy (on-resonance) and 9.6 fb _1 approximately 40 MeV below 
this energy (off-resonance). 

Charged-particle momenta are measured in a tracking system consisting a five-layer double- 
sided silicon vertex tracker (SVT) and a 40-layer central drift chamber (DCH), both immersed in a 

7 Charge conjugation is implied here and throughout this paper unless explicitly stated. 



8 



1.5 T axial magnetic field. BABAR achieves an impact parameter resolution of about 40 /jm for the 
high momentum charged particles in the B decay final states, allowing the precise determination 
of decay vertices. The tracking system covers 92% of the solid angle in the CM frame. 

Charged-particle identification is provided by measurements of the energy loss {dE/dx) in the 
tracking devices (SVT and DCH) and by an internally reflecting ring-imaging Cherenkov detector 
(DIRC) covering the central region. A K-ir separation of better than four standard deviations (a) is 
achieved for momenta below 3 GeV/c, decreasing to 2.5 a at the highest momenta in the B decay 
final states. Photons are detected by a CsI(Tl) electromagnetic calorimeter (EMC). The EMC 
provides good energy and angular resolution for detection of photons in the range from 20 MeV to 
4 GeV. The energy and angular resolutions are 3% and 4 mrad, respectively, for a 1 GeV photon. 

3 EVENT SELECTION 

Hadronic events are selected based on track multiplicity and event topology. We fully reconstruct 
B° — * p + p~ candidates from their decay products p^ — > ir^ir and 7r° — * 77. Charged track 
candidates are required to originate from the interaction point, and to have at least 12 DCH hits 
and a minimum transverse momentum of 0.1 GeV/c. Charged pion tracks are distinguished from 
kaon and proton tracks via a likelihood ratio that includes dE/dx information from the SVT and 
DCH for momenta below 0.7 GeV/c and the DIRC Cherenkov angle and number of photons for 
higher momenta. They are also distinguished from electrons primarily on the basis of the EMC 
shower energy and lateral shower moments. 

We reconstruct ir° mesons from pairs of photons, each with a minimum energy of 30 MeV, with 
the shower shape consistent with the photon hypothesis, and not matched to a track. The typical 
width of the reconstructed ir° mass is 7 MeV. We accept ir° candidates in the invariant mass interval 
±15 MeV from the nominal mass. We select p candidates to satisfy 0.52 GeV < m n7T < 1.02 GeV. 
The helicity angles 9\ and 62 of p + and p~ are defined as the angles between the 7r° direction in the 
p rest frame and the direction of the p boost with respect to the B as shown in Fig. |2I We restrict 
the helicity angles to the region —0.75 < cos < 0.95 to suppress combinatorial background and 
reduce acceptance uncertainties due to low- momentum pion reconstruction. 




Figure 2: Definition of helicity angles (9± and 62) in B° — ► p + p decays. 



To reject the dominant continuum background (from e + e — * qq events, q = u,d,s,c), we 
require | cos 9t\ < 0.8, where 9t is the angle between the thrust axis of the B candidate and the 
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thrust axis of the rest of the tracks and neutral clusters in the event, calculated in the CM frame. 
The distribution of | cos 6t\ is sharply peaked near 1 for jet-like events originating from qq pairs and 
nearly uniform for the isotropic decays of the B meson. We also construct a Fisher discriminant 
(J-) that combines 11 variables: the polar angle of the B momentum vector, the polar angle of the 
-B-candidate thrust axis, both calculated with respect to the beam axis in the CM frame, and the 
scalar sum of the CM momenta of charged particles and photons (excluding particles from the B 
candidate) entering nine coaxial angular intervals of 10° around the i?-candidate thrust axis |1U| . 

We identify B meson candidates using two nearly independent kinematic variables J2J, the beam 
energy-substituted mass toes = [(s/2 + Pi • pb) 2 /E^ — p#] 1//2 and the energy difference AE = 
{EiEb — Pi • Pb — s/2)/i/s, where (E^pj) is the e + e~ initial state four-momentum, and (Eb,Pb) 
is the four-momentum of the reconstructed B candidate, all defined in the laboratory frame. For 
signal events the toes distribution peaks at the B mass and the AE distribution peaks near zero. 
Our initial selection requires toes > 5.2 GeV and \AE\ < 0.2 GeV, while the signal resolution is 
roughly 3 MeV and 50 MeV, respectively. The selected sample contains 54042 events most of which 
are retained in sidebands of the variables for later fitting. 

4 ANALYSIS METHOD 

We use an unbinned, extended maximum-likelihood (ML) fit to extract simultaneously the signal 
yield and angular polarization. There are three event yield (nj) categories j: signal, continuum qq, 
and BB combinatorial background. We define the likelihood for each B° — > p + p~ event candidate i: 

3 

C i = ^2n j V j (x i ;d), (1) 
i=i 

where each of the Vj{xi,a) is the probability density function (PDF) for measured variables x,. 
The numbers a parameterize the expected PDFs of measured variables in each category. We allow 
for multiple candidates in a given event by assigning to each selected candidate a weight of 1/iVj , 
where iVj is the number of candidates in that same event. The average number of candidates per 
event is 1.27. The extended likelihood for a sample of N can ^ candidates is 

(3 \ JVcand /ln£\ 
~H n i i n ex p ■ ( 2 ) 

The seven fit input variables x% are toes, AE, T , invariant masses of the p + and p~~ candidates, 
and the corresponding helicity angles 6\ and 02- The correlations among the fit input variables for 
the data and signal Monte Carlo (MC) JJ| are found to be small (typically less than 5%), except 
for angular correlations in the signal as discussed below. The PDF, Vj{xi \ a), for a given candidate 
i is the product of those for each of the variables and a joint PDF for the helicity angles, which 
accounts for the angular correlations in the signal and for detector acceptance effects. We integrate 
over the angle 4> between the two decay planes shown in Fig. leaving a PDF that depends only 
on the two helicity angles and the unknown longitudinal polarization fraction /x = Tl/T, where 

and r are the longitudinal and total decay widths. The differential decay width j2j is defined 

as: 

¥a a =7(7(1-/L)sin 2 g 1 sin 2 g 2 + / L cos^ 1 cos 2 2 j . (3) 

r d cos 6*i d cos 62 4 L 4 J 
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The PDF parameters a, except for f^, are extracted from MC simulation and on-resonance 
rriES an d AE sidebands. The MC resolutions are adjusted by comparisons of data and simulation 
in abundant calibration channels with similar kinematics and topology, such as B — > Dtt, Dp 
with D — > Kirir, Kir. To describe the signal distributions, we use Gaussian functions for the 
parameterization of the PDFs for tbes an d AE, and a relativistic P-wave Breit-Wigner distribution 
for the p resonance masses. The angular acceptance effects are parameterized with empirical 
polynomial functions for each helicity angle and are included in the joint helicity-angle PDF as a 
product with ideal distribution in Eq. Q. For the background PDFs, we use polynomials or, in the 
case of rriES; an empirical phase-space function ^21- I n the background PDF we incorporate a small 
linear correlation between the curvature of the phase-space function and the event shape variable T . 
The background parameterizations for the p candidate masses also include a resonant component 
to account for p production in the continuum. The background helicity-angle distribution is also 
separated into contributions from combinatorial background and from real p mesons, both described 
by polynomials. The PDF for T is represented by a Gaussian distribution with different widths 
above and below the mean for both signal and background. 

There is a fraction of incorrectly reconstructed (fake) B° — ► p + p~ events expected in the 
selected sample. This happens when at least one candidate photon in a tt° candidate or one 
charged track in a p candidate belongs to the decay tree of the other B. MC simulation shows 
that about 30% of selected B° — > p + p~ events with longitudinal polarization do not have the 
correct decay tree reconstructed, while about 20% of the events have both correctly and incorrectly 
reconstructed decay candidates. We do not account explicitly for the fake events in the signal PDF 
parameterization since they introduce substantial tails in the distributions; these tails have weak 
discrimination power from background and we cannot fully rely on the MC simulation of the fake 
component distributions. We account for these effects in the signal efficiency evaluation. 

MC simulation indicates that about 5% of the events in the final sample are from other B 
decays. This background, arising mainly from b — * c transitions, is explicitly accounted for in the 
fit. PDFs for this background are taken from MC including a contribution from charmless decays 
such as B —>■ pit, p°p + , pK*, a±7r, and a\p. The branching fractions for these and many other 
exclusive modes were taken from the most recent experimental measurements |13| or extrapolated 
from other results with flavor symmetry approximation. Their contribution was shown to be well 
accounted for by a single ^-background fit component. In this analysis we do not explicitly include 
a fit component for other partial waves, including non-resonant decays, with the same final-state 
particles selected within the p resonance mass window. These types of decays are assumed to be 
negligible, they are significantly suppressed by the mass and helicity-angle information in the fit, 
and they are examined in the mass and helicity-angle distributions as discussed below. 

The event yields rij and polarization are obtained by minimizing the quantity x 2 = — 21n£. 
The dependence of x 2 on a fit parameter rij or fi is obtained with the other fit parameters floating. 
We quote statistical errors corresponding to a unit increase in x 2 ■ The statistical significance of 
the signal is defined as the square root of the change in x 2 when constraining the number of signal 
events to zero in the likelihood fit. 

5 PHYSICS RESULTS 

The results of our maximum likelihood fits are summarized in Tabled The statistical significance 
of the B° — > p + p~ signal is 5.9 a. We find that the decay amplitude is predominantly longitudinal. 
To compute the branching fraction, we assume equal production rates for B°B° and B + B~ . To 
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check the stability of our results we refit removing each variable from the fit and find consistent 
results. The number of fitted events, statistical significance, branching fraction and polarization 
measurement errors, and the ML fit x 2 value are well reproduced with generated MC samples. 

The projections of the fit input variables are shown in Fig. 13 The projections are made after 
a requirement on the signal-to-background probability ratio V s \ g /Vhkgi where V s ig and Pbkg are 
the signal and background PDFs defined in Eq. l|Tj). but with the PDF for the plotted variable 
excluded. The histograms show the data with about 40—60% of signal retained, the lines show 
the PDF projections from the full sample. Both mass and helicity-angle projections provide no 
evidence for non-resonant B decays with the same four-pion final states. 

To check the sensitivity of our results to the presence of non-resonant B° — > p^ir^-ir and 
B° — > ir + iT~ir ir decays, we explicitly include a fit component assuming phase-space decay model. 
The selection requirements alone suppress the B — > pirir (47r) efficiency by one (two) order (s) of 
magnitude relative to B° — > p + p~ ■ The fit results with non-resonant component are consistent 
with our assumption of negligible non-resonant contribution. We exclude the hypothesis that all 
the signal is non-resonant B — > pmr (47r) with 4.7 a (5.4 a) statistical significance. However, we 
cannot exclude that B — > pirn signal could represent (9 ± 15)% (statistical errors only) of our 
nominal B° — > p + p~ event yield, where we ignore interference effects. 

Table 1: A summary of the fit results. The efficiency includes systematic errors, the significance 
with systematic uncertainties is quoted in parentheses, while for other results, the first error is 
statistical and the second systematic. 



Reconstruction efficiency 


3.9^| % 


Signal event yield (n s j g ) 


93j|? ± 9 


Statistical significance 


5.9 g (5.3 cr) 


Branching fraction (B) 


(27 ±1 t 5 7 ) x 10- 6 


Signal polarization (f£) 


0.99 ±g;{£ ± 0.03 



6 SYSTEMATIC STUDIES 

Systematic uncertainties in the ML fit originate from assumptions about the PDF parameters. 
Uncertainties in the PDF parameters arise from the limited statistics in the background sideband 
data and signal control samples. We vary them within their respective uncertainties, and derive 
the associated systematic error on the event yield (9%). The signal remains statistically significant 
under these variations (5.3 cr including systematics). Additional systematic errors in the number of 
signal events originate from the uncertainty in the other S-decay cross- feed (3%) that was studied 
with exclusive MC generated samples. 

The systematic errors in the efficiency are due to track finding (2% for two tracks), particle 
identification (2% for two tracks), and 7r° reconstruction (13% for two 7r°s). The efficiency in the 
ML fit to signal samples, calculated as the ratio of the fit signal yield over the number of fully 
reconstructed decays in the fit sample, is less than 100% because of fake combinations passing 
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Figure 3: Projections onto the variables m^s, AE, m^+^o, m^-^o, cos 6\, and cos 62 after a require- 
ment on the signal-to-background probability ratio 'P s i g /'Pbkg with the PDF for the plotted variable 
excluded. The histograms show the data, the solid (dashed) line shows the signal-plus-background 
(background only) PDF projection. 
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the selection criteria. We account for this in the efficiency evaluation and assign a systematic 
uncertainty of 7% taken to be 1/2 of the inefficiency. The reconstruction efficiency depends on 
the decay polarization. We calculate the efficiencies using the measured polarization and assign a 
systematic error ( + },2%) corresponding to the total polarization measurement error. Smaller sys- 
tematic uncertainties arise from event-selection criteria, MC statistics, and the number of produced 
B mesons. 

For the polarization measurement, we also include systematic errors from PDF variations that 
account for uncertainties in the detector acceptance and background parameterizations (0.025). 
The biases from the resolution in helicity-angle measurement and dilution due to the presence of 
the fake combinations are studied with MC simulation and are accounted for with a systematic 
error of 0.02. 

7 SUMMARY AND DISCUSSION 

We have observed the decay B° — » p + p~ , measured its branching fraction B = (27l 6 ly) x 10 -6 and 
determined the longitudinal polarization fraction fi = 0.991q'q7±0.03. Our results are preliminary. 
This completes the measurements of the isospin-related B — > pp modes El- The measured 
branching fraction agrees well with the more recent predictions in the range of (18-35) x 10~ 6 4 , 
and the suppression of the transverse amplitude (by a factor of m p /mB) was expected (!]. These 
measurements improve our understanding of the dynamics of hadronic weak decays and allow 
experimental tests of effective theories and factorization (2JII3II]- 

The rates of the B° — > p + p~ and B + — > p°p + decays appear to be larger than the corresponding 
rates of B — > -kit decays ^3], while the recent measurement of the B + — > p°K* + branching frac- 
tion |Hj does not show significant enhancement with respect to B — * irK decays |13j . This indicates 
that the penguin pollution in the B —* pp decay is smaller than in the tttt case, as predicted prior 
to the experimental measurements of the B — > pp modes |Hj- The dominance of the CP-even lon- 
gitudinal polarization in both B° — > p + p~ and B + — ► p°p + decays will also simplify CP-violation 
studies. 

At the same time, a relatively smaller experimental limit on the decay rate can be achieved in 
the B° — > p°p° decays than in the B° — ► decays, this allows better experimental limit on the 
amount of penguin pollution in the B — > pp decay amplitudes. Since the tree contribution for the 
B° — > p°p° decay is color-suppressed, the decay rate is sensitive to the penguin diagram in Fig. ^ 
Using the earlier BABAR measurements ,6.. we obtain a 90% confidence level upper limit on the 
ratio of the longitudinal amplitudes Al in the B — > pp decays: 

\A L (B° -> pV)| 2 + \A L (B° p°p°)\ 2 _ B(B ^p p )xf L (B ^p°p ) ^ n in m 
2 x \A L (B+ -^p°p+)\ 2 

In the above calculation we assume conservatively the P u — ► p [) p^ decay polarization to be fully 
longitudinal {fi = 1), use the average branching fraction measurements for the B and B decays, and 
assume \Al(B + — ► p°p + )\ = \Al(B~ — ► p°p~)\ with only a tree-diagram contribution. A similar 
experimental limit for the B — > and B + — > 7r + 7r° decay amplitudes is 0.61 |14j . though 

still restrictive on the penguin pollution in the measurement of a 5]. The above observations 
make the B° — > p + p~ decay a promising channel to study CP-violation and set constraints on the 
weak-interaction angle a. 



B(B- 



14 



8 ACKNOWLEDGMENTS 



We are grateful for the extraordinary contributions of our PEP-II colleagues in achieving the 
excellent luminosity and machine conditions that have made this work possible. The success of 
this project also relies critically on the expertise and dedication of the computing organizations 
that support BABAR. The collaborating institutions wish to thank SLAC for its support and the 
kind hospitality extended to them. This work is supported by the US Department of Energy and 
National Science Foundation, the Natural Sciences and Engineering Research Council (Canada), 
Institute of High Energy Physics (China), the Commissariat a l'Energie Atomique and Institut 
National de Physique Nucleaire et de Physique des Particules (France), the Bundesministerium fur 
Bildung und Forschung and Deutsche Forschungsgemeinschaft (Germany), the Istituto Nazionale 
di Fisica Nucleare (Italy) , the Foundation for Fundamental Research on Matter (The Netherlands) , 
the Research Council of Norway, the Ministry of Science and Technology of the Russian Federation, 
and the Particle Physics and Astronomy Research Council (United Kingdom). Individuals have 
received support from the A. P. Sloan Foundation, the Research Corporation, and the Alexander 
von Humboldt Foundation. 

References 

[11 M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, 652 (1973). 
|21 G. Kramer and W.F. Palmer, Phys. Rev. D 45, 193 (1992). 
R. Aleksan et al, Phys. Lett. B 356, 95 (1995). 



[3 
[4 

[5 

[6 

[7 

[8 
[9 
[10 

[11 

[12 
[13 
[14 



D. Ebert, R.N. Faustov, and V.O. Galkin, Phys. Rev. D 56, 312 (1997); A. Ali, G. Kramer, 
and C.-D. Lu, Phys. Rev. D 58, 094009 (1998); Y.-H. Chen et al, Phys. Rev. D 60, 094014 
(1999); H.-Y. Cheng and K.-C. Yang, Phys. Lett. B 511, 40 (2001). 

M. Gronau and D. London, Phys. Rev. Lett. 65, 3381 (1990); Y. Grossman and H. Quinn, 
Phys. Rev. D 58, 017504 (1998). 

BABAR Collaboration, B. Aubert et al, BABAR-PUB-03-018, |hep-ex/0307026| submitted to 
Phys. Rev. Lett. 



BELLE Collaboration, J. Zhang et al, BELLE-2003-6, hep-ex/0306007 ( submitted to Phys. 
Rev. Lett. 

BABAR Collaboration, B. Aubert et al, Nucl. Instrum. Methods A479, 1 (2002). 

PEP-II Conceptual Design Report, SLAC-R-418 (1993). 

CLEO Collaboration, D.M. Asner et al, Phys. Rev. D 53, 1039 (1996). 

The BABAR detector Monte Carlo simulation is based on GEANT: R. Brun et al, CERN 
DD/EE/84-1. 

ARGUS Collaboration, H. Albrecht et al, Phys. Lett. B 241, 278 (1990). 
Particle Data Group, K. Hagiwara et al, Phys. Rev. D 66, 010001 (2002). 
BABAR Collaboration, B. Aubert et al, Phys. Rev. Lett. 91, 021801 (2003). 



15 



